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ABSTRACT 

We explore isothermal shock formation in non-equatorial, adiabatic accre- 
tion flows onto a rotating black hole, with possible application to some active 
galactic nuclei (AGNs). The isothermal shock jump conditions as well as the 
regularity condition, previously developed for one-dimensional (ID) flows in the 
equatorial plane, are extended to two-dimensional (2D), non-equatorial flows, 
to explore possible geometrical effects. The basic hydrodynamic equations with 
these conditions are self-consistently solved in the context of general relativity 
to explore the formation of stable isothermal shocks. We find that strong shocks 
are formed in various locations above the equatorial plane, especially around a 
rapidly-rotating black hole with the prograde flows (rather than a Schwarzschild 
black hole). The retrograde flows are generally found to develop weaker shocks. 
The energy dissipation across the shock in the hot non-equatorial flows above 
the cooler accretion disk may offer an attractive illuminating source for the re- 
processed features, such as the iron fluorescence lines, which are often observed 
in some AGNs. 

Subject headings: transonic accretion, isothermal shock waves - black hole physics, 
adiabatic flows, hydrodynamics - relativity 



1. Introduction 

A realistic model of the AGN central engine may include the effects of the magnetic field. 
The black hole magnetosphere is studied first by Blandford & Znajek (1977) in the context 
of winds and jets from radio-loud AGNs. The work has been extended to magnetospheric 
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physics of accreting AGNs by various authors (e.g., Phinney 1983; Takahashi et al. 1990; 
Tomimatsu & Takahashi 2001; Takahashi et al. 2002, hereafter TRFT02; Rilett et al.2004, 
in preparation). In this case, plasma particles should be frozen-in to the magnetic field 
lines, and hence the accreting fluid should fall onto the black hole from regions above the 
equatorial plane along the field lines(see, e.g., Figure 2 of TRFT02; Figure 1 of Tomimatsu 
& Takahashi 2001). Then, magnetohydrodynamics (MHD) should become important to 
describe the motion of the particles associated with the background field. The resulting 
relativistic MHD shocks are explored by TRFT02 and Rilett et al. (2004, in preparation). 

In general, the MHD shocks can be hydro-dominated or magneto-dominated (Takahashi 
2000, 2002, hereafter T02). Obviously the MHD case is, however, very complicated. There- 
fore, the main motivation of our current paper is, as a starting point, to investigate the 
hydrodynamic limit, which should be valid in the case of small magnetization. Here, we 
adopt a model which should apply to the hydro-dominated shocks where the magnetic field 
does not make a significant contribution to the properties of the shocks under a weak field 
limit, because in such a case hydrodynamics should primarily control the shock formation. 

One-dimensional (ID), hot accretion flows around a black hole, generally treated as 
ideal hydrodynamical fluid, have been investigated by various authors (Sponholz & Molteni 
1994; Kato et al. 1996; Kato, Fukue, & Mineshige 1998). It has been found that such 
accretion flows must be transonic, and will become supersonic before reaching the event 
horizon, while it is subsonic at infinity. Once such a fluid becomes supersonic, it is likely 
that a standing shock wave will develop when shock conditions are met. There are roughly 
three types of standing shocks: adiabatic (Rankine-Hugoniot) shocks, isothermal shocks, 
and isentropic compression waves (Abramowicz & Chakrabarti 1990). The first attempt was 
made by Yang & Kafatos (1995) to self-consistently study the relativistic isothermal shock 
formation around black holes. In the case of isothermal shocks, the postshock fluid can lose 
substantial energy and entropy across the shock while the fluid temperature is continuous 
across the shock location (Lu & Yuan 1997). Lu & Yuan (1998, hereafter LY98) examined 
the isothermal shock formation in one-dimensional (ID) adiabatic hot flows in the Kerr 
geometry, for various flow parameters including black hole rotation. 

In our current paper, we extend the work by LY98 on isothermal shock formation in ID 
adiabatic flows in the equatorial plane, to two-dimensional (2D) calculations for flows above 
the equatorial plane, to investigate geometrical effects. One of our major motivations is to 
explore the possibility that shocks produced in such flows act as a high energy radiation 
source for some reprocessed features, such as the iron fluorescent lines, which are observed 
from some AGNs. It is generally considered that a hot illuminating source above the cooler 
disk is required to produce such features (Fabian et al. 1989). 
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Section §2 introduces our basic equations and assumptions. The hydrodynamic fluid 
equations are solved in the context of general relativity We discuss the isothermal shock 
conditions and the stability of shocks. The results are presented in Section §3 where we 
display shocks for various representative values of angular position, fluid energy, angular 
momentum, and black hole spin. Discussion and concluding remarks are given in the last 
section, §4. 



2. Basic Equations & Assumptions 

2.1. Background Kerr Geometry 

The space-time is assumed to be stationary (d/dt = 0) and axially-symmetric (<9/<90 = 
0) around a rotating black hole. The background Kerr metric is then expressed in the 
Boyer-Lindquist coordinates as 

, 2 / 2Mr\ , 2 4aMrsin 2 e , , , 
ds 2 = - ( 1 — J dt 2 dtdcf) 

A sin 2 9 , ,2 X , 9 ^ 
+ - d(f) 2 + -dr 2 + Zd6 2 , (1) 

where A = r 2 — 2Mr + a 2 , £ = r 2 + a 2 cos 2 9, A= (r 2 + a 2 ) 2 — a 2 A sin 2 9 with the metric sig- 
nature being (— , +, +, +). Throughout this paper, the distance is normalized by the gravita- 
tional radius r g = GM/c 2 where c, G, M being the speed of light, gravitational constant, and 
the black hole mass, respectively. The hole's event horizon is then = (1 + Vl — a2 ) r g- No- 
tice that 9 = 9 s h = 7r/2 corresponds to the shocks in the equatorial plane while 9 = 9 s h ^ n/2 
denotes the non-equatorial shocks discussed here. 



2.2. Hot, Adiabatic Fluid above the Equatorial Plane 

Chakrabarti (1996a,b), Lu et al. (1997) and LY98 have studied general relativistic, 
equatorial flows in the Kerr geometry, whereas many of the previous works were based on the 
Pseudo- Newtonian potential (Paczynski & Wiita 1980), which is not capable of reproducing 
the frame-dragging effect. To examine the two-dimensional (2D), relativistic hot accretion 
flow, here we assume an ideal, Boltzmann gas in a non-equatorial plane. Such an accreting 
fluid spirals around the symmetry axis due to angular momentum and accretes onto the black 
hole. The poloidal path of the fluid is assumed to be conical (u e = 0) since the preshock fluid 
slowly starts falling from a distant location as a subsonic flow. That is, the fluid has only 
the radial velocity (u r < 0) and the azimuthal velocity (u^ ^ 0). Therefore, flow particles 
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will spiral around the black hole rotation axis onto the hole. The schematic diagram of such 
an accreting fluid is displayed in Figure 1. The non-equatorial accreting fluid is conically 
falling onto the hole (i.e., #=constant and u r < in the poloidal plane), spiralling around 
the symmetry axis in the azimuthal plane. 

In this model perturbing forces to the fluid that might cause any acceleration in the 
6 direction is assumed to be negligible, keeping the flow at constant 6 in a conical flow. 
This assumption is justified because we are exploring the weak field limit of MHD flows, 
which should apply to flows with small magnetization. As noted in Section §1 already, fluid 
particles in such flows should be frozen-in to the field lines and flow along the filed lines. 
Therefore, the flow geometry in our case is governed by the magnetospheric structure, rather 
than by the hydrostatic equilibrium of non-magnetic fluid. In this sense it is not appropriate, 
in our case, to adopt the approach of conventional, non-magnetic 2D thick accretion disk (or 
torus) models. Especially, it may be emphasized that our version of the 'conical model' is 
quite different from the conical equilibrium flow often adopted in conventional non-magnetic 
thick disk models. 

Solving the geometry of the magnetospheric structure in a realistic manner is an ex- 
tremely difficult problem, and only a very simple approximate version has been carried out 
(see, e.g., Tomimatsu and Takahashi 2001). Therefore, in our current paper which should ap- 
ply to the weak field limit, we use, as our first approximation, the flow geometry we adopted 
already in our previous MHD accretions model where the infalling plasma is frozen-in to 
the radial field lines in a conical geometry, and hence the flows in different 9 directions are 
decoupled from each other (see TRFT02, Rilett et al. 2004, in preparation). 

Following the earlier works on relativistic accretion shocks, particularly LY98, we assume 
that a dynamical time-scale of the accretion process is much shorter than that for the energy 
(or thermal) dissipation during the fluid accretion. The fluid obeys the equation of state for 
an ideal gas 

P = -^pT , (2) 

fJjpUlp 

where the primary thermodynamic property of the fluid is characterized by the locally mea- 
sured temperature T and the thermal pressure P. p, p p and m p are the rest-mass density of 
the fluid, the mean molecular weight of the composite particles and the mass of a particle, 
respectively. k B is the Boltzmann constant. The polytropic form is adopted as 

P = Kp~< = Kp 1+1/n , (3) 



where the adiabatic index 7 is constant, and the polytropic index n is correspondingly 
constant, too. We use 7 = 4/3 (or n = 3) for our relativistic fluid. K is a measurement of 
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the entropy S of the gas where S = c v logK (c v being a specific volume heat), expressed also 
via equations (2) and (3) as 

K = — — -T . (4) 

We consider the fluid to be stationary and axisymmetric. The fluid trajectory is conical in 
the poloidal plane and is spiralling around the polar axis onto the hole. Thus, #=constant 
or u e = 0, but u r < and vfi ^ 0. Due to the space-time symmetry, specific total energy E 
and axial angular momentum component L of the fluid are conserved along the shock-free 
fluid's path and are written as 

E = -pu t , (5) 
L = fiu,p , (6) 

where \i = (P+e) / p is the relativistic enthalpy of the fluid, and e = p+nP is the total energy- 
density including the internal energy term. With the use of the four-velocity normalization 
u ■ u = — 1 where u = (V, u r , 0, u^) is the four- velocity of the fluid, we find 

l + u r u r + (u t ) 2 V eff (r,£) = , (7) 

where V e ff(r,£) = g tt — 2£g t ^ + £ 2 g^ is called the effective potential (Lu, Yu, & Young 
1995) with g a @ being the Kerr metric tensor components. £ = Lj E = —u^/ut is the specific 
angular momentum of the fluid which is conserved along the whole flow if viscous dissipation 
is negligibly small (weak viscosity limit). Hence, we get u t in terms of V e ff(r,£) as 

1 + u r u r ^ 



On the other hand, the definition of the local sound speed is given by c s = (dP/de) 1 ^ 2 
following LY98, which is then rewritten as 



^4=(i + ;)^(S, (») 



de \ n J \de 
and by definition of the enthalpy we also have 



p = ^^ = l + (l + n)Kp 1 / n . (10) 
P 



Thus, combining equations (9) and (10), we find 



9 f dP\ 11 

\-nc\ = \-n — =- ——, = - , 11 

V de J 1 + (1 + n)Kp l l n p : v 1 
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where we now express the local sound speed in terms of the enthalpy (or vice versa). Ac- 
cordingly, therefore, equation (5) can be rewritten as 

*=(i±gfV*«). (12) 

It is useful to notice here that both E and L will change across the shock location in the 
case of isothermal shock, but they change in such a way that the ratio L/E (i.e., £) will stay 
unchanged (see a subsequent section). Notice that K remains the same along the shock-free 
fluid due to the adiabatic assumption (no heat in/out), but does change across the shock. 
There is another constant along the fluid's path called mass accretion rate M. The spherical 
mass accretion rate is 4nr 2 (—u r )p in the radial direction, and thus the mass accretion rate 
per unit spherical area by incoming fluid is given by 

M = r 2 (-u r )p , (13) 

which is constant along the whole flow. Since we are interested in shock formation where 
entropy is generated, let us follow LY98 and introduce another conserved quantity called the 
entropy accretion rate 

2 1 11 

n ° s 2 , (-« r ), (14) 

7 (1 -ncf)_ 



M = K n M = r 



where M is conserved in a shock-free fluid, but changes in a shock-included fluid due to the 
entropy generation (i.e., K) at the shock location. 



2.3. Regularity Conditions 

It is known that any black hole accretion must be transonic (Chakrabarti 1990a, 1996a), 
so it is important to locate the critical points (or sonic points to a particular observer) to 
make sure that the fluid is physically acceptable. Taking a derivative of equation (12) with 
respect to r, we obtain 

du r N 

-¥ = d> (15) 



where 



N ^ (u r ) 2 dg rr 2c 2 s l_dVefi (iQ) 

2(1 +u r u r ) dr r 2V eff dr ' 1 ' 



C 2 U r 



D = — — . (17) 

u r 1 + u r u r v ' 



-7- 



In equation (15), D(N = 0) = is required (regularity conditions) for du r /dr to be finite. 
Therefore, at critical points we get 



= — | (is) 



d{lng rr ) ' 
-dr 4/r 



C 5 



(19) 



where the subscript "c" denotes the critical points. Notice that the right-hand-side of equa- 
tion (18) is the radial three- velocity component of the fluid measured by a corotating observer 
(this will be explained later). In other words, the critical points are equivalent to the sonic 
points only for such an observer. Multiple critical points (i.e., an inner point r* n , an inter- 
mediate point r™ d and an outer point r° ut ) are known to exist in general for a set of (E, £) 
(e.g., Chakrabarti 1990b). This means that the location of the critical points will be different 
before and after the shock in a shock-included fluid since E changes across the shock. The 
fluid must pass through a critical point on both sides of the shock. It is worth while to note 
that in isothermal shocks the postshock fluid cannot pass through the same critical points 
for the preshock fluid due to the energy dissipation at the shock location. This will be ad- 
dressed more in a later section. Depending on the critical point that the fluid goes through, 
kinematic and thermodynamic profiles of the fluid will be different. In the framework of our 
model, however, it is not terribly important through which critical point the accreting fluid 
will pass as far as it produces a dissipative shock near the black hole. Therefore, we will 
only investigate the preshock fluid that passes through the outer critical point r° ut and the 
postshock fluid that passes through the inner critical point r l c 



c ' 



2.4. Isothermal Shock Conditions 

In our scenario the adiabatic fluid forms a standing shock (i.e., stationary shock location) 
somewhere between the preshock outer critical point r° ut and the postshock inner critical 
point r*™ on its way to the horizon. According to Abramowicz & Chakrabarti (1990), stand- 
ing shocks in general can be categorized into three types: (l)adiabatic shocks or Rankine- 
Hugoniot shocks, (2)isentropic compression waves and (3)isothermal shocks. In the case of 
(1), the fluid by definition does not release any energy across the shock, carrying the gen- 
erated entropy and its thermal energy with it (Fukue 1987; Chakrabarti 1990b; Lu et al. 
1997). Thus, radiative cooling mechanism is extremely inefficient with the thermal energy 
being advected. In type (2), the shock radiates an energy equivalent to the generated entropy 
at the shock such that the entropy remains unchanged at the shock (Chakrabarti 1989). In 
the case of type (3), a fraction of the preshock fluid's energy and the entropy generated 
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are lost from the fluid's surface at the shock location such that the temperature (therefore 
the sound speed) remains continuous across the shock (Lu & Yuan 1997, LY98). In other 
words, the dissipative cooling processes are so efficient that the energy is not advected with 
the fluid. A realistic shock could be between these two extremes (adiabatic and isothermal). 
However, in this paper we investigate case (3), because it is very interesting in the sense that 
a large amount of energy releases can be expected in the vicinity of the central black hole, 
which may be related to the observed emission features in AGNs. 



Since the fluid total energy decreases across the isothermal shock, we have 

1/2 



/ 1 + U{u rl \ , 2 v 



E 2 
E s h 



( 1 + u r 2 u r2 \ 1/2 . , 

V -Veff J /[ 



nc 



sh 



) . 



Ei — E 2 , 



(20) 

(21) 
(22) 



where the subscript "1" and "2" denote a preshock and postshock quantity at the shock 
location, respectively. c s h = c s (r sh ) is the local sound speed at the shock location. From the 
entropy accretion rate, we have 



Mi 



sh 



sh 



■'sh 



7(1 - nc 2 sh ) 



"sh 

M l - nc lh). 

Due to the shock formation, the fluid becomes hotter and entropy is generated at the shock 
front (i.e., mathematical discontinuity with zero-thickness). Thus, M. increases instanta- 
neously at the shock location. However, the radiative cooling process may be so efficient 
that a fraction of the total entropy may be lost as radiation from the fluid surface. There- 
fore, the postshock temperature remains continuous and the postshock entropy becomes 
smaller than that of the preshock. That is, T\ = T 2 but Ai± > Ai 2 although Mi = M 2 . The 
radial momentum flux density is given by 



-u\ 



-u r 2 ) 



(23) 
(24) 



+ u r 



From conservation of momentum flux density in the radial direction, it is found that 



7«rl 



^sh 
lU r2 



(25) 



(26) 



Based on the above jump conditions in addition to the regularity conditions, there are 9 
equations (12), (14), (18), (19), (20), (21), (23), (24) and (26) for 9 unknowns r* n , r° ut , r sh , 



- 9- 



E 2) Mi, M 2) c s hi u\ and u 2 . Once a set of (a, 9 s h, E\, £) is specified, then the corresponding 
shock location r s h and all the other unknown quantities are uniquely determined. For the 
sake of completeness, it should be noted that any physically acceptable shocks in a transonic 
accretion flow must satisfy the boundary conditions (r* n < r sh < r° ut ). 

From equations (20) and (21), we note that E\ > E 2 at the shock location r s h because 
| till > \u 2 1. That allows some energy dissipation E sh . Equations (23) and (24) also imply 
that the entropy accretion rate should decrease across the shock for the same reason. Now, 
we may extract some interesting qualitative implications from these equations. For example, 
we note that more energy release E s h should naturally be expected from strong shocks 
because a relatively large jump in the radial velocity u r across such shocks is made possible. 
Secondly, we clearly see the general trend for the effect of the direction, 9 sh , on specific 
angular momentum of the fluid, £. For instance, £ cannot be large near the polar region 
since the centrifugal barrier is correspondingly small compared with that near the equator. 
Then, the non-equatorial fluid will have smaller rotational velocity u^. Consequently, the 
flow could gain its radial component u r more efficiently than its azimuthal component u^, 
which allows the preshock flow to become supersonic more quickly (or sooner) in the radial 
direction. The resulting shock location, therefore, will be further from the black hole for 
small 9 s h as compared with the shocks near the equator (large 9 s h) ■ Furthermore, the energy 
dissipation E s h from such (more distant) shocks will be small, due to the fact that only a 
small amount of fluid's gravitational potential energy can be released, as opposed to a large 
amount of energy conversion from the gravitational potential in regions close to the hole. 
We will come back to this issue in Section §3. 

As mentioned in the earlier section, the critical points r c are determined by the fluid 
parameters (E, £). In dissipative shock-included flows, the preshock fluid energy E x (and M.\) 
decreases across the shock, which will subsequently determine the postshock critical points 
corresponding to different E. In other words, the postshock fluid with (E 2 ,Ai 2 ) cannot 
pass through the critical points already determined by the preshock flow with (Ei,Aii). 
Instead, new critical points for the postshock flow will be determined by the postshock flow 
parameter (E 2 , Ai 2 ). Therefore, one must find these new critical points based on (E 2 , Ai 2 ) for 
the postshock flows to describe global shock-included transonic flows. Such a flow topology 
will be illustrated in §3.5. 

2.5. Stability Analysis 

Although possible shock locations are found from the jump conditions, they need to be 
stable against a small perturbation. Otherwise, the shock will disappear as a result of a small 
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change in its location. To check the stability, we perturb the radial momentum flux density 
(equivalent to pressure) by a small change in the shock location r s h, a method originally 
developed by Yang & Kafatos (1995). Thus, a small perturbation of the momentum flux 
density is 



When the momentum flux density is not balanced on both sides of the shock, the shock is 
not stable. Suppose that a shock location is slightly perturbed by Sr. If A > when 5r > 
(or Sr < 0), the postshock momentum flux density becomes larger (or smaller) for accretion. 
Consequently, the shock location will be further shifted towards increasing (or decreasing) r, 
and the fluid will never find a stable shock location. If A < 0, on the other hand, the shock 
will always go back to the original location satisfying momentum equilibrium regardless of 
the sign of Sr. Therefore, A < is required for stable shocks for accretion. 



It is interesting to see how flow quantities are changing as the accreting gas falls onto 
the black hole. Here, we will estimate some of the important flow quantities for global, 
shock-included, transonic fluids as a function of radial distance r. First, we explore the 
kinematics of the flows through the four- velocity components u r , and the angular velocity 
Q = u^/u 1 as well. The physical flow velocity can be described by the radial three- velocity 
v r = [u r u r /(l + u r u r )] 1 / 2 in a corotating reference frame (CRF) which is a locally-flat space. 
v r is known to be equal to the speed of light at the horizon. The Mach number is defined in 
CRF by M = v r /c s . The relativistic enthalpy of the flow can be written as fi = 1/(1 — nc 2 s ) 
in the unit of fluid rest-mass energy (see equation (11)). The thermal pressure P can be 
derived from equation (3) as P oc p 7 . Then, the local temperature of the flow goes as 
T oc P/p using equation (2). By the mass conservation law, we get the local density of the 
flow as p oc l/(r 2 \u r \). Here, we define a normalized pressure P = P/K = p 7 , temperature 
T = P/p = k B T I '(fj, p m p ) and the mass density p = p/M = l/(r 2 |-u r |), respectively. 




(27) 



2.6. Characteristics of Global, Shock-Included Flows 



3. 



Results 



In this section we present our results and consider their physical implications. The 
obtained isothermal shock solutions are displayed for various parameters. Based on previous 
AGN observations and theoretical speculations (e.g., Iwasawa et al. 1996a,b; Wilms et al. 
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2001), the mass of the central supermassive black hole is chosen to be M = 1O 7 M , and the 
hole is in general rapidly-rotating (i.e., an extreme Kerr black hole). We choose a = 0.99 
in most cases. However, in Section §3.4 the black hole spin dependence is investigated by 
adopting various values of a. The accretion fluid is rotating with the Keplerian angular 
velocity of Qk = l/(a + r 3 / 2 ). We consider both direct flows (ai > 0) and retrograde flows 
(ai < 0). 

As seen by LY98, we find stable, multiple shock locations (r*^, where r l Jl < r^f ) for 
the same set of parameters (a, 9 s h, Ei, i) satisfying all the requirements (i.e., shock jump con- 
ditions, boundary conditions, and stability analysis). The preshock flow which does not reach 
the event horizon is often called a-type flow. In order for a-type preshock flows to become 
physically acceptable, there must exist a shock. However, it is still uncertain through which 
shock location, between r™ h or r^f , the preshock flow will "transit" to the postshock flow 
(i.e., degeneracy of shocks). There may be some other physical conditions (for instance, vis- 
cous dissipation or magnetic field) that uniquely determine which shock is actually realized, 
but this is beyond the scope of this paper. To avoid complexity caused by the degeneracy, we 
will include only the inner range of shock locations r*^ (which is closer to the hole). Although 
we are aware of the existence of the "degeneracy", throughout this paper the outer range 
of shock locations are not considered. The reason is that, as mentioned already, our work 
is originally motivated by the X-ray observations of the "reprocessed features" such as the 
iron lines from some AGNs. These lines are considered to come from a cool accretion disk 
illuminated by a hot high energy radiation source above the disk plane. The whole system 
should, however, be very close to the central black hole(Fabian et al. 1989), because these 
lines often exhibit broadened features with a large gravitational redshift(Nandra & Pounds 
1994). Therefore, we are primarily interested in the inner shock formation close to the hole. 



3.1. Isothermal Shock Locations above the Equatorial Plane (0° < 9 sh < 90°) 

Solving the basic equations (conservation laws + shock conditions) introduced earlier, we 
found various stable shock locations for different fluid parameters such as (9 sh , Ei,i). The fol- 
lowing four angles 9 s h are considered: 5° (near-axis region), 30° (lower polar region), 60° (upper 
equatorial region) and 80° (near equatorial region). The energy-dependence of the shock lo- 
cations for various E\ is also examined in the following manner: E 1 = 1.005, 1.01, 1.015 and 
1.02 for comparison where E\ is expressed in the unit of the particle's rest-mass energy m p c 2 . 
To avoid complications, we fix the black hole spin at a = 0.99 (with the prograde flows). 
First, the cross-sectional distribution in the (X,Y)-plane of the obtained shock locations is 
displayed in Figure 2 where E 1 is (a)1.005, (b)1.01, (c)1.015 and (d)1.02 for the four selected 



- 12 - 



9 sh : 5°, 30°, 60° and 80° from top to bottom except for (d). As explained earlier, the 
outer range of shock locations is not shown there. It should be reminded that a different 
shock location r s h is accompanied by a different angular momentum £ in this figure. We will 
discuss this later in §3.3. 

We find that the non-equatorial shock develops further away as £ increases for a fixed 
Ei and 9 s h, consistent with LY98 in the equatorial case. Both the inner critical point r*™ 
and the outer critical point r° ut will shift inward as £ increases for a fixed E\ and 9 s h- There 
exist a minimum £ (£ m in) and a maximum £ (£ m ax) for each set of (Ei,6 s h) for which the 
physically valid, stable shock formation is possible. In other words, the shock is unable to 
develop because of the violation of the boundary conditions and/or because of the instability 
under ceratin £ with (E 1: 9 sh ). This is also commonly found in the equatorial shocks. The 
relation between the shock location r s h and the angular momentum £ for a = 0.99 is seen in 
Figure 3 where r sh is plotted against £ for various angle 9 sh . E\ = 1.005, 1.01, 1.015 and 
1.02 from top to bottom except for 9 sh = 5° where E x = 1.005, 1.01 and 1.015. Notice that 
infinitely many shocks do continuously exist between any two consecutive shocks plotted 
here although only representative shocks are shown (in dots) in this figure. 

The range of the non-equatorial shock locations Ar s h in general tends to become nar- 
rower as E\ increases for a fixed 8 s h, where Ar s h = r™£ x — r^ n . r^ m and r^ ix are the 
innermost (smallest) shock location and the outermost (largest) shock location in the inner 
branch, respectively. From Figure 2, such a trend is obvious. For a fixed angle 9 s h, the 
average location of the inner critical point f* n = {r™ min + r t l ! l max )/2 seems to be larger while 
the average location of the outer point f° ut = [r°^ nin + r°^ ax )/2 seems to be smaller as E\ 
increases, where the additional subscript "min" and "max" denote the minimum value and 
the maximum value of the critical points for a fixed E 1 and 9 s h, respectively. Another inter- 
esting feature with regard to the angular dependence is that Ar sh roughly becomes larger 
with decreasing 6 sh for a fixed E 1 . That is, both the innermost shock location r™ n and the 
outermost shock location r™£ x occur farther away from the hole when the accreting flow is 
nearer the polar region although r™£ x shifts outwards more than r™ ra . This is because the 
average r* n and the average r° ut both tend to become smaller with increasing 9 s h- However, 
the degree of decrease in r° ut is much larger than that of r* n , which forces A s h to be smaller 
with larger 9 sh . Interestingly, no shocks can form when E 1 = 1.02 with 9 sh = 5° as seen 
in (d). This is due to the presence of the outer critical points very close to the black hole. 
With such a small r° ut (together with a certain value of r* n ), there can be only a narrow 
(radial) region between r* n and r° ut . In the case of 5° in (d), this region becomes so narrow 
that no jump conditions can be met within this range, resulting in no shock formations at 
all regardless of the angular momentum £. 
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As discussed in the earlier section §2.4, the average shock location seems to shift out- 
wards with decreasing 9 sh , which is qualitatively predicted as a result of a quick gain of 
the radial velocity due to smaller angular momentum £. In general, therefore, if the total 
preshock energy E\ is relatively large when 9 s h is small, the formation of shocks may not 
necessarily be expected. In other words, the formation of shocks in a polar region may 
require less energetic hydrodynamic accretion fluid. The characteristic transition of these 
important locations are summarized in Table 1. 

3.2. Properties of Dissipative Shocks 

Figure 4 shows a three-dimensional plot where the shock location r s h, the shocked fluid 
angle 9 s h and the shock strength M1/M2 are displayed for various E-y. We take a = 0.99 
(for the prograde flows) with E\ = 1.005, 1.01, 1.015 and 1.02 for (a), (b), (c) and (d), 
respectively. Here, Mi/M 2 is the fluid Mach number ratio just before the shock and just 
after the shock. For clarity purpose, the range of the plot varies from (a) to (d). Once again, 
no shock is found for 9 sh = 5° when E\ = 1.02 (see panel (d)). For a global geometry of the 
whole shock location, refer to Figure 2. 

One of the easily noticeable features in this figure is that Mi/M 2 increases with decreas- 
ing r s h down to a certain distance (peak radius r p ) and then turns to decrease inside r p for 
a fixed Ei. This is again the same pattern found in LY98 for 9 = 90° case. Such a trend 
can be seen more explicitly as 9 s h increases (compare 5° and 80°, for instance). However, as 
Ei increases, the rate at which Mi/M 2 increases outside r p seems to become smaller while 
the rate at which M 1 /M 2 decreases inside r p appears to be larger for a fixed 9 s h- That is, 
the slope of the Mi/M 2 — r s h curve tends to become smaller with increasing E\ (compare 
(a) 1.005 and (d)1.02, for instance). Therefore, we can conclude that the preshock fluid with 
a small (large) E\ and a large (small) 9 s h is more likely to produce a strong (weak) shock 
for a fixed shock location r sh . 

The energy dissipation E^/Ei (or the ratio of the energy released to the preshock 
energy) from the preshock fluid is shown (in percent) as a function of r s h and 9 s h in Figures 5. 
The panels (a)-(d) correspond to the same energy E\ as in Figure 4. This ratio basically 
follows the same pattern as the Mi/M 2 except for a specific detail. Since the E sh jE x — r sh 
curve qualitatively resembles the Mi/M 2 — r sh curve, the energy dissipation E sh becomes 
larger as the shock occurs closer to the hole until the peak radius r p if it exists at all. The 
release of the energy from the hot flow tends to be more efficient when 9 s h is large for a fixed 
energy E\, and the variation of E s h/E\ can be nearly an order of magnitude if the shock is 
developed near the disk plane (i.e., 9 sh ~ 90°) which can be seen in (a) and (b) in Figure 5. 
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This can be qualitatively interpreted in terms of the energy conversion of the fluid, as we 
predicted in the section §2.4. The Shock developed further away can dissipate only a small 
gravitational potential energy whereas the shock close to the hole can convert more potential 
energy for dissipation. The maximum energy release can be as high as ~ 18% whereas the 
minimum energy release seems less than 1% in the case of a = 0.99 (with prograde flows). 
As another trend, more energy E s h tends to be released from a hot flow with a smaller E\. 
Overall, it can be concluded that more energy release is expected from a low-latitude (i.e., 
large 9 sh ) hot flow with a small energy E\. 



3.3. Angular Dependence 

In order to see the angular dependence of the specific angular momentum of the shocked 
fluid, £, Figure 6 shows the £ versus 9 s h relation. We choose preshock fluid energy and black 
hole spin at fixed values, E 1 = 1.005 and a = 0.99, respectively, in order to better understand 
the effects due to 9 sh alone. We find that at each 9 sh , £ ranges from £ min to £ max for different 
radial shock location r s h, (although it may be difficult to see this in the figure because 
£ m x — £min is much smaller than the variation of £ over 9 s h) ■ One can express the angular 
momentum as a function of r — sh, as £ = £{r s h). Note that LY98 also finds similar behavior 
(i.e., the r s h dependence of £). In fact, this is a general behavior found for any other shock- 
included global solutions. On the other hand, we note that for shock-free global accreting 
flows the r dependence disappears, and hence a single value of £ should be determined by 
specifying 9 alone. 

We clearly see (in Figure 6) the general trend that the fluid angular momentum decreases 
with decreasing 9 s h- That is physically consistent with the fact that the centrifugal force (due 
to the gravitational potential well) decreases towards the polar axis. Thus, the fluid must 
possess smaller angular momentum in the polar region (i.e., for small 9 s h) for accretion to be 
realized. A similar trend has been found also in the work of the MHD shocks by TRFT02. It 
may be possible, in principle, to analytically show the dependence of the physically allowed 
angular momentum £ as a function of 9 s h, for a given value of r s h (i.e., £(r s h, 9 s h))- However, 
that is beyond the scope of our present work. For our computation purposes these two 
variables (£(r s h) and 9 s h) are, to start with, treated as independent of each other. However, 
the valid solutions numerically found do indicate that £ is uniquely determined for a given 
set of the shock location (r sh ,9 sh ). 
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3.4. Black Hole Spin Dependance 

Here we explore the effect of black hole spin on shock properties. To decouple the 
effect of the hole spin from the other effects, we hold the fluid energy and the angle at 
fixed values of E 1 = 1.01 and 9 sh = 30°, respectively. Figure 7 demonstrates the energy 
dissipation E sh /E 1 (in %) as a function of the spin parameter a and shock location r sh . We 
take a = —0.99, —0.5, 0, 0.5 and 0.99. Since the fluid angular momentum I is always taken 
to be positive here, the cases a = —0.99 and —0.5 refer to retrograde flows. 

First we note that energy is most efficiently dissipated around a rapidly-rotating black 
hole in the presence of prograde flows. The retrograde flows with a = —0.99, on the other 
hand, appear to produce the least energy release, which corresponds to small Mach number 
ratio M 1 /M 2 . Due to the frame- dragging by the hole, the prograde flows can quickly acquire 
more velocity (especially the rotational component) in the course of its accretion, as opposed 
to the retrograde case where the fluid rotation appears even smaller (in the local inertial 
frame). Therefore, the prograde flows can afford more change in the kinetic motion at a 
shock location, and thus it may be "easier" for the prograde flows to dissipate more energy 
using the obtained kinetic motion. This viewpoint can be directly translated into the large 
Mach number ratio M 1 /M 2 that we generally observe in the case of the prograde accretion. 
This may naturally explain large energy release E sh jE x at shocks that we normally find. On 
the contrary, the retrograde flows generally cannot afford as much change as the prograde 
can in terms of the kinetic motion, and therefore only weak shocks can develop with small 
Mach number ratio (or small energy dissipation). 

Keeping in mind that we are only interested in the inner range of the shocks here, 
the minimum shock location r™ n tends to shift radially outwards with decreasing spin a. 
The maximum shock location r™£ x is the greatest when a = 0. This suggests that in the 
retrograde flows no shock formations should be expected in a region very close to the hole, 
as opposed to the prograde case. 

We can conclude from Figure 7 that very strong shocks (with more energy release) are 
expected in the inner region for the prograde flows. 

3.5. Global Shock-Included Accreting Flow Solutions 

As explained in Section §2.3, in the case of isothermal shocks, the global shock-included 
transonic solution must pass through the outer critical (sonic) point determined by the 
preshock parameters and simultaneously goes through the inner critical (sonic) point de- 
termined by the postshock parameters. Before looking into the detailed of particular flow 
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dynamics, let us illustrate such a flow topology for a specific example. Figure 8 shows the 
Mach number M of the transonic flow as a function of the radial distance r (in logarithmic 
scale) for a shock-included transonic solution with a = 0.99, £ = 1.463 and 9 s h = 30°. The 
preshock flow with (Ei, = (1.01,5.27 x 10 -5 ) goes through the outer critical point at 
r out _ 53 2r g , then makes the transition via the shock at r s h = llr 5 into the postshock flow 
with (E 2 , M2) = (0.98258, 1.06 x 10 -5 ) passing through the inner critical point at r*™ = 2.58r ff 
before reaching the horizon. Notice that the preshock flow topology is completely indepen- 
dent of the postshock flow topology (both drawn in dark curves) in that they individually 
possess their own critical points (i.e., non-mutual points), and the flow transition from the 
former to the latter is allowed only through the shock (denoted in arrow). The other flows 
(shown in light gray curves) with different values of M. are unphysical solutions in terms of 
the fact that the jump conditions are not satisfied and that no critical points exist. 

Here, we present three types of global, shock-included, transonic, accreting flow solutions 
in the case of a = 0.99 (with prograde flows): Figure 9 for flow 1, Figure 10 for flow 2 and 
Figure 11 for flow 3. The flow parameters are tabulated in Table 2 for each flow. Each figure 
displays the radial component \u r \ and the azimuthal component of the four velocity, 
the angular velocity Q (solid curve) with the equatorial Keplerian velocity Qkep (dotted 
curve), the radial three- velocity in CRF |t> r | (solid curve) with the local sound velocity c s 
(dotted curve), Mach number M, specific enthalpy p, flow temperature T, local density p and 
thermal pressure P respectively. The enthalpy here is also normalized by the fluid rest-mass 
energy m p c 2 . The outer/inner critical (sonic) point is denoted by a large dot. The value 
at the event horizon is shown by another large dot. Either a vertical arrow or a small dot 
shows the transition from a preshock flow to a postshock flow through a stable isothermal 
shock. Note that rhor < r T < r sh < r° ut for physically acceptable shocks. In Table 3, some 
representative fluid quantities for each flow is displayed in physical units. To make concrete 
calculations, we adopt the following relation between the (dimensionless) mass accretion 
rate for the non-equatorial hot gas ifn hot ) and the observed net accretion rate ifn net ) in the 
units of the Eddington mass accretion rate (~ 1.4 x 10 25 g/s) for a 1O 7 M black hole mass. 
Here, we take rhhot ~ 0.01m net and m net = 0.1 (relevant for radio-quiet AGNs) based on the 
assumption that the most of the accretion rate (99%) is contributed from the cool accreting 
gas in the accretion disk, whereas only a fraction of the net accretion rate (~ 1%) is due to 
the hot non-equatorial accreting fluid. 

As seen in the earlier sections, flow 1 in Figure 9 goes through a relatively strong shock 
(M 1 /M 2 ~ 5.4) as a consequence of the dramatic increase in the density p and the pressure P. 
The flow temperature T seems to be monotonically rising up to a peak point r p where the local 
sound speed is the maximum. In fact, all the thermodynamic quantities, such as p, T and 
p, are all correlated to the sound speed. Notice that the local sound speed c s is continuous 
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across the isothermal shocks because the flow temperature T is by definition continuous. 
Passing through r p , the sound velocity starts decreasing which triggers a relatively large 
decrease in the thermodynamic properties such as p, p, P and T. The angular velocity 
Q = u^/u 1 roughly follows the Keplerian value (in the equator) almost all the way down 
to the horizon. The radial three- velocity |t> r | appears to become maximum just before the 
shock and then transits through a strong shock, rising up again rapidly to become supersonic 
before entering the horizon. It is seen that, as expected, |t> r | = 1 at the horizon. 

Flow 2 in Figure 10 shows somewhat different behaviors compared to flow 1. First, the 
shock strength is weaker (M1/M2 ~ 3.0). In terms of the kinematics, the angular velocity 
reaches a maximum value at some point and then slows down in the azimuthal direction 
from this point on. Notice that Q(r) > for its entire trajectory. An interesting difference 
here is that the radial velocity is still decreasing after the shock for about ~ 5r g and then 
turns to pick up the (radial) speed. This is probably due to a relatively large shock location 
(r s h ~ 13r 3 ). Thermodynamic quantities also behave in different ways here. After the shock, 
the flow temperature does not drop in this case. Instead, it continues to rise gradually just 
after the shock and the rate of the change becomes larger towards the horizon. This behavior 
is originating from the sound speed c s (r). The rest of the thermodynamic quantities (p, p 
and P) also follow the same pattern. 

The shock strength for flow 3 in Figure 11 is somewhere between flow 1 and flow 2 
(Mi/M 2 ~ 4.8). An outstanding feature in this case is a large deviation of Q(r) from the 
Keplerian angular velocity in the equator. Since the shock location is relatively close to the 
hole, the postshock (radial) speed increases immediately after the shock. In both temperature 
and density profiles, a strange behavior can be seen after the shock. In temperature profile, 
the postshock fluid temperature appears almost constant producing a plateau-like pattern, 
whereas density profile p shows a step-like change in p. Again, such a feature of these 
thermodynamic quantities can be explained by the sound speed c s (r). 

It appears in general from the above three cases that the postshock fluid does not become 
supersonic right away if the shock is formed at a relatively distant location (say, r sh > 10r g ), 
in which case the radial velocity can be still slowing down. It turns out, however, that it is 
not easy to classify all the possible shock-included flows according to the angle 9 s h because 
the flow energy E\ as well as the angular momentum I are all certainly related to the flow 
dynamics, which allows a varieties of flow dynamics even for a fixed 9 s h- For this reason, we 
will not try to simply attribute the above hydrodynamic/thermodynamic features in each 
flow to just one parameter out of (Ei,£, 9 s h). 
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4. Discussion & Concluding Remarks 

Although our current work is partially motivated by the work of LY98, it also originates, 
in important ways, from our very recent work on black hole magnetospheres in accretion- 
powered AGNs - TRFT02 and Rilett et al. (2004, in preparation), where we explored 
adiabatic, relativistic MHD shocks produced in the black hole magnetosphere. We showed 
that strong shocks can indeed be formed in such flows for various relevant choices of flow pa- 
rameters. In the current paper we extended these previous studies to explore the relativistic 
hydrodynamic flows which should apply to the case of weak magnetization. The reason is 
that the physics involved in the exact MHD case would be far more complicated, as noted 
already by TRFT02. Owing to these complexities, it was not straightforward for us to study 
the exact global shock-included MHD accretion flow solutions further in detail in a wider 
parameter/solution space. 

As our next step, in the current paper we assumed a simple model of conical accretion 
flows, due to the fact that modelling more realistic flows (such as the ones in hydrostatic 
equilibrium) would be very complicated, especially in the framework of the relativistic non- 
equatorial accretions considered here. An appropriate force balance under the general rel- 
ativistic geometry should be taken into account in more sophisticated models, but that is 
beyond the scope of our present work. We, however, would like to stress that our results still 
represents (at least qualitatively) important physical characteristics of shock formation in 
non-equatorial accretion flows. Also, in the presence of the magnetosphere, the fluid particles 
would be frozen-in to the field and hence flow along the field lines. For such a situation, 
application of conventional thick accretion disk (or torus) models is not appropriate. 

We considered only the inner range of the shock formations in regions relatively close 
to the central engine, because of our interests in some X-ray observations of the reprocessed 
emission, such as the iron fluorescence lines, from some AGNs. Our results may offer the 
possibility of a high energy source in various broad regions (0° < 6 s h < 90°) above the disk 
plane. However, we find that the strongest shocks (high Mi/M 2 with large E s h) should 
develop near the equator (large 6 s h), although the quasi-polar shock (small 6 sh ) is also 
possible. We find no shocks in the polar region (9 sh = 5°) when the preshock fluid energy is 
relatively large although the shock-free accretion is physically allowed. The magnitude of the 
energy release from the shock roughly increases as the shock location gets closer towards the 
black hole, as already found by LY98 by their ID studies. In addition, however, we further 
find, from our 2D calculations, that the shock-induced energy release E s h greatly depends, 
not only on the fluid energy E\ and angular momentum £, but also on the angle of the shock 
location 9 s h- The average shock strength (thus energy release) tends to be weaker towards 
the polar region. 
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Although TRFT02 considered adiabatic shocks and hence no energy dissipation, in the 
current paper we adopted the isothermal shocks because that ensures a substantial amount 
of energy release at the shock locations. Also, through the stability analysis our energy 
source (i.e., the shock) is found to be stable. Although we adopted adiabatic flows in our 
current work, it may be noted that Das, Pendharkar, & Mitra (2003) recently explored 
the isothermal shock formations in the isothermal flows, by adopting pseudo-Schwarzschild 
gravitational potentials and ID equatorial flows. These authors concluded that their shocks 
also will release substantial amount of energy, which could be physically sufficient to become 
a radiation source for a strong X-ray flare. 

In the present work, it is argued that a single isothermal shock formation between the 
two sets of accreting flows (preshock/postshock flows) is very likely with a substantial amount 
of energy dissipation. One may also ask whether it is possible to have a sequence of shock 
formations one after another in the course of the accretion. For example, a preshock flow with 
energy E 1 gets shocked at a shock location r sh ^ becoming a subsonic postshock flow with 
energy E2. The same flow with energy E2 then becomes supersonic and develops another 
shock at r s h y 2 where r s h^ < fsh,i, making the transition to another subsonic postshock flow 
with energy E 3 where E 3 < E 2 < E 1 . Such a sequential shock formation (namely "shock 
cascade" or "shock avalanche" ) could be a very interesting phenomenon. Although discussion 
of such a possibility is beyond the scope of the present work, it is interesting, in a future 
work, to investigate such a possibility in relation to some observational applicability. 

Before closing, we emphasize that a major justification for the choice of our version 
of a conical flow geometry follows from our recent work of TRFT02 on magnetospheric 
MHD accretion flows. However, here we adopted the hydrodynamic flows, as a limiting case 
of weak magnetization, in order to take advantage of the current fortunate situation that 
the relativistic hydrodynamic (equatorial) shocks in ID accreting flows have already been 
extensively studied by many authors (e.g., Chakrabarti 1989; Abramowicz & Chakrabarti 
1990; Lu et al. 1997; Lu & Yuan 1997, and LY98). Our current investigation, however, is 
new and valuable in the sense that we have explored 2D non- equatorial shocks in a fully 
relativistic manner, with a possible application, for instance, to offer an attractive definite 
source for the energy dissipation in regions very close to the black hole. 

We are grateful to Darrell Rilett and Masaaki Takahashi for enlightening suggestions 
on our model. KF also thanks Maki Fukumura for computational assistance. ST thanks 
colleagues in Institute of Astronomy, Cambridge, especially Drs. A. Fabian, G. Miniutti, 
M.J. Rees, and K. Iwasawa, for valuable discussions and comments. We are especially 
indebted to the anonymous referee for providing a number of constructive suggestions to 
improve the manuscript. 
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Rotational Axis Non-Equatorial 

A Accreting Hot Fluid 




BH Accretion Disk 



Fig. 1. — A side view of the schematic geometry in our model. A non-equatorial accreting 
fluid is conically accreting in the poloidal plane (u 6 = 0) onto the black hole, spiralling 
around the rotational axis (u r <0,u^ ^ 0) with a constant angle 9 s h- It passes through the 
outer critical point r° ut , gets shocked at r s h and passes through the inner critical point r* n 
before reaching the horizon. The shock location r s h increases with increasing specific angular 
momentum of the fluid £. See the text for details. 
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Fig. 2. — The range of isothermal, stable shock locations around a rapidly-rotating black 
hole (a = 0.99 with prograde flows) for various angle 8 s h for fixed energy E\. The black 
hole's event horizon is = 1.141r 5 (central dot). No shock is found for E\ = 1.02 with 
s h = 5° in (d). Note that only the inner branch of shock locations are shown. 




Fig. 3. — The shock location r s h versus angular momentum £ for various E\ and 6 s h in the 
case of a = 0.99 (with prograde flows). The dots represent the obtained shock locations. 
Each curve corresponds to the respective curves in Figure 2. From top to bottom, E\ = 
1.005, 1.01, 1.015 and 1.02 except for 6 s h = 5° where no shock is found for E\ = 1.02. 
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Fig. 4. — Shock strength (Mach number ratio) M\jMi as a function of the shock location 
r s h and the angle 9 s h for the same shocks as in Figure 2 (r™ n < r s /j < r^ ix ). The dots 
represent the obtained shocks. E = 1.005, 1.01, 1.015 and 1.02 for (a), (b), (c) and (d), 
respectively. We take a = 0.99 with prograde flows. 
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Fig. 5. — The ratio of the energy dissipation to the preshock energy E s h/E\ (in %) for the 
same shocks as in Figure 2. 
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Fig. 6. — The specific angular momentum of the flow I versus the angle 9 s h- We take 
8 s h = 5°, 30°, 60° and 80°. Each dot is a collection of many smaller dots that correspond to 
particular angular momentum at different radial shock locations r s h- The other parameters 
are a = 0.99 and E 1 = 1.005. 



-28- 




Fig. 7. — The energy dissipation ratio E s h/Ei (in %) as a function of the shock location r s h 
and the black hole spin a for E\ = 1.01 and 9 s h — 30°. The dots represent the obtained 
shocks. We choose a = —0.99, — 0.5, 0, 0.5 and 0.99. Negative black hole spin cases are 
considered to have retrograde flows (i.e., a£ < 0). 
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Fig. 8. — A flow topology of a global shock-included accreting flow with E\ = 1.01 and 
6 s h = 30° for a = 0.99. A preshock transonic flow with (Ei, M.\) passing through an outer 
critical point r° ut develops an isothermal shock at r s h and becomes a subsonic postshock flow 
with (E2, M.2)- The flow then passes through an inner critical point r* n becoming supersonic 
again before reaching the horizon. The physically valid preshock/postshock flows through 
the critical points are denoted by dark curves, while light grey curves refer to unphysical 
flows. 
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Fig. 9. — Global shock solutions for an adiabatic fluid (flow f) with E\ = 1.02, 1 = 
2.104, = 2.005r 9 for 9 s h = 80°. The horizon, inner and outer critical (sonic) points 
are respectively denoted by large dots from left to right in each figure. The shock is repre- 
sented by either a small dot or a vertical arrow. From the upper-left panel to the upper-right 
panel, the radial component \u r \ of the flow, the azimuthal component of the flow and 
the angular velocity Q (solid curve) with the Keplerian value Qkep [dotted curve) are shown. 
From the middle-left panel to the middle-right panel, the radial three- velocity |t> r | of the 
flow in the corotating reference frame (CRF) (solid curve) with the local sound velocity c s 
(dotted curve), the Mach number M and the relativistic enthalpy fi are displayed. From the 
lower-left panel to the lower-right panel, the fluid local temperature T, the density p, and 
the thermal pressure of the flow P are plotted. For details, see Tables 2 and 3. 
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Fig. 10. — Global shock solutions for an adiabatic fluid (flow 2) with E\ = 1.015, 1 = 
1.4325, r s h = 12.75r s for 6 sh = 30°. The notations are the same as in Figure 9. 
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Fig. 11. — Global shock solutions for an adiabatic fluid (flow 3) with E% = 1.005, 1 = 
0.2958, r s h = 6.95r 9 for 9 S h = 5°. The notations are the same as in Figure 9. 
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Table 1. Some Important Locations as a Function of E 1 and 6 s h- 
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1.005 


80° 


3.625 


18.00 


14.38 


1.262 


146.0 


1.01 


5° 


6.188 


20.32 


14.13 


4.238 


57.71 


1.01 


60° 


4.050 


12.85 


8.810 


1.420 


69.28 


1.01 


80° 


3.330 


10.80 


7.471 


1.278 


71.05 


1.02 


30° 


5.062 


10.62 


5.559 


3.048 


25.62 


1.02 


60° 


2.060 


7.800 


5.741 


1.515 


31.77 


1.02 


80° 


1.650 


7.000 


5.352 


1.316 


33.56 



Note. — The length (or distance) is in the unit of r g . See 
text for notations, a = 0.99 for all cases. 
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Table 2. Global Shock-Included, Transonic, Non-Equatorial Accretions for Various Flow 

Parameters. 



Flow 


Ei 


I 


dsh 


r sh/r g 


r in /rg 


r out /rg 


M 1 /M 2 




-M 2 (10- 5 ) 


Figure 


1 


1.02 


2.104 


80° 


2.005 


1.347 


33.65 


5.398 


15.43 


2.551 


9 


2 


1.015 


1.4325 


30° 


12.75 


2.697 


37.85 


2.987 


9.265 


3.046 


10 


3 


1.005 


0.2958 


5° 


6.95 


4.059 


133.3 


4.761 


1.906 


0.3896 


11 



Note. — a = 0.99 for all cases. The black hole mass M is 10 7 M Q . 
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Table 3. Thermodynamic Properties of Accreting Hot Flows at the Shock Location. 



Flow 


Ei 




r sh/r g 


c s /c 


ni(10 10 ) a 


n 2 (10 10 ) a 


Pi(10 5 ) b 


P 2 (10 5 ) b 


Figure 


1 


1.02 


80° 


2.005 


0.2420 


12.06 


73.00 


79.67 


482.0 


9 


2 


1.015 


30° 


12.75 


0.1351 


0.4051 


1.232 


0.8335 


2.535 


10 


3 


1.005 


5° 


6.95 


0.1258 


1.253 


6.129 


2.235 


10.93 


11 



a in 1/cm 3 
b in dyne/cm 2 

Note. — The subscripts "1" and "2" denote the preshock and postshock quantities, respectively. We 
adopt black hole mass M = 10 7 Mq, and m^ot = 0.01rrt ne t = 10~ 3 . Sec the text for notations. 



